Predictions of Quark-Gluon String Model for pp at LHC 
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Abstract 

Soft multiparticle production processes in hadronic col- 
lisions are considered in the framework of the Quark- 
Gluon Strings Model and the model predictions are com- 
pared with data from SppS and Tevatron. Predictions for 
LHC energies are given. 



1 Introduction 

In hadronic collisions soft interactions are dominant. 
These are characterized by large distances (r ^ 1/Aqcd) 
and the pQCD cannot be used for their study. 

Regge-pole theory is the main method for description 
of high-energy soft processes and the Pomeron is the 
main object of this approach, which is associated with 
cylinder-type diagrams of expansion in QCD. 
expansion is a dynamical expansion and the speed of 
convergence depends on the kinematical region of study- 
ing process. At very high energies many terms of the 
expansion (multi-pomeron exchanges) should be taken 
into account. 

The models of Quark-Gluon Strings (QGSM) T and 
Dual-Parton Model (DPM) 2 , are based on nonpertur- 
bative notions, combining expansion in QCD with 
Regge theory and using parton structure of hadrons. 
Here we concentrate on QGSM, which has been suc- 
cessfully applied to many different problems of strong 
interactions: hadronic mass spectrum [3], widths of res- 
onances [4j, relations between the total cross-sections, 
residues of Regge poles [4], behavior of hadronic form 
factors [5], baryon number transfer [S], and multiparti- 
cle production at high energies [H [H [71 [8] . 

In this paper we briefly describe QGSM, compare its 
predictions with SppS and Tevatron data on charged par- 
ticles pseudorapidity and multiplicity distributions and 
give predictions for LHC. 



2 QGSM 

In QGSM the production of a particle is defined through 
production of showers and each shower corresponds to 
the cut-pomeron pole contribution in the elastic scatter- 
ing amplitude. For hadron-hadron collisions the cross- 
section of n cut-pomereon exchange (among arbitrary 
number of uncut pomerons) is calculated in "quasi- 
cikonal" approximation and has the following form [9]: 



n — 1 / 

Z 



1 - exp{-z} ^ 



1=0 



II 



(1) 



Where z = C-//Xexp{A^}, \ = + a'p^, and f = 
ln(s/so). The values of parameters 7, A, and a'p 
which characterize the residue and the trajectory of the 
pomeron are found from fit to data on pp and pp total in- 
teraction and elastic scattering cross-section in Ref. [TU] : 
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2.14 GeV-2, i?2 ^ 3 3 GeV-2, A = 0.12, a'p 



0.22 



GeV~^. The parameter C, which is related to small- 
mass diffraction dissociation of incoming hadrons during 
the rescattering is equal to 1.5 [llj. The sum of a„ {n 
= 1,2,3, . . .) defines the cross-section of non-diffractive 
interaction. 

The association of the pomeron with cylinder-type di- 
agrams leads to the fact that in a single cut-pomeron 
diagram there are two chains of particles (strings). Anal- 
ogously, in case of n cut-pomerons there are 2n chains. 
Inclusive distribution in rapidities and multiplicity dis- 
tribution of particles for non-diffractive events can be ex- 
pressed in terms of rapidity distributions, /n(C,J/), and 
multiplicity distribution, Wn{£,,N), for 2ri-chains [T]: 



dy ^ 

a(C,Af)=^a„(OW^„(C,iV)- 



(2) 



In order to calculate the distribution of hadrons pro- 
duced during the fragmentation of the strings one should 
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take into account that the hadron can be produced in 
each of the 2n-chains and the (di-)quarks , which stretch 
the strings, carry only a fraction of energy of the incom- 
ing protons. In these terms, the inclusive spectra can 
be written as convolutions of the probabilities to find a 
string with certain rapidity length and the fragmenta- 
tion functions, which defines the distribution of hadrons 
in the string breaking process. For pp collisions the func- 
tion fn{£,,y) is written as follows [I]: 



+2{n^l)F^}-J{x+)F'j::\x^) 



where x± — {l/2)[^/x^^+x^ ± x], xt — 2m'^/^ and 
X is the Feynman-a; of produced hadron h. is the 
density of hadrons h produced at mid-rapidity in a sin- 
gle chain and its value is determined from experimental 
data. In fact these are the only free parameters in the 
model that are fixed from fit to data. In articles [l],[6] 
and |7] the following values are found for them: = 
0.44, = 0.055, qP = 0.07. The first two terms in 
Eq. (I3|) correspond to the chains, which connect "va- 
lence" quarks and di-quarkf0, and the last term corre- 
sponds to chains connected to the sea quarks-antiquarks. 
Here we do not consider the contribution of sea di-quarks 
and anti-diquarks, their contribution is important in the 
spectra of anti-baryons produced in the forward region 
[12]. The functions F^^"\x) being a convolution of the 
structure ip{x), and fragmentation G{x) functions are 
writen as follows: 



^.1"H^±) = 3 / dxyju.,^,{n,x)Gt (^) 

1 X 

+ -J dxV'd„„,(n,x)G^' (^) , 
J^^J"H^±) = IJ^ dx^ud{n,x)Gt, (^) 
+ ^ / dx^/juu{n, x)G'^^^ (^^) , 
F,^i:H^±) = [ dx^^_{n,x)Gt (^) 



'In case of pp collisions they must be replaced by 



Where S is the strangeness suppression parameter (S ~ 
1/3). 

In the model, the structure and fragmentation func- 
tions are determined by the corresponding Regge asymp- 
totic behaviors in the regions x — s- and a; — s- 1 and for 
the full range of x an interpolation is done. 

The structure functions of a proton are parameterized 
)as follows [ilZ]: 

' (|]„(a;,n) = =Gda:-"«(l-x)"«-2"«+", 

^ud{x, n) = G„d:r"--2a« „ 

^uuix,n) = GuuX (1 — x) , 

Mx,n) = C^a;-"«(i _a;)2("«-"B)-"*.+"-i. 

Where aii=0.5, Q!b=-0.5. Ci are determined from nor- 
malization condition: 



order to have white objects required by QCD. 



-,h{n) 



h{n)i 



This change is in 



'ipiix, n)dx = 1. 



General technique of constructing fragmentation func- 
tions is presented in Ref. 13J. For instance, for frag- 
menting to a pion the parameterizations are (see [1] and 



Gf{z) = Gl{z) = {l^zf-^^-+^ 
Gl\z) = Gl{z)^{l^z)-^-+^ 
Gf{z) = Gr{z) = {l-zf--+^ 
Gtd (^) = Gl, (z) - (1 - z)"-+^-2"« (1 _ z + ^V2) 
GlUz)^{l-zr-+^-^"- 



We do not list fragmentation functions for kaons and 
(anti-)protons, they are taken from Ref. 017]. 

In hadronic interactions diffractive processes play an 
important role. Densities of produced particles in these 
processes are small, especially in the central rapidity re- 
gion, but their cross-section is not negligible with re- 
spect to non-diffractive cross-section and they must be 
taken into account in calculations of characteristics of 
secondary particles. In Ref. [14] single- and double- 
diffractive processes are described in terms of dressed 
triple-reggeon and loop diagrams. For instance, for cal- 
culating the rapidity distribution of particles produced 
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in soft single-difFractive events the first equation in ^ 
must be replaced by (see [8]): 



da 
dy 



E 



(4) 



Here da/dC, is the derivative of single-diffraction dissoci- 
ation cross-section, C, — M'^ / s, and M is the mass of the 
diffracted system. 

The model does not contain "odderon"-type singular- 
ities with negative signature, which could lead to a dif- 
ference of pp and pp scattering, so at energies > 100 
GeV all characteristics oipp and pp interactions coincide. 

3 Numerical results 

QGSM has been used successfully to describe data on 
secondary hadron inclusive cross-sections integrated over 
transverse momentum, such as rapidity and multiplicity 
distributions, ratios of particles and it explains many 
characteristic features of hadron-hadron soft interac- 
tions. But for calculating pseudorapidity distribution 
from Eq. ([2]) or ([4]) we must know mean transverse mo- 
mentum of particles. 

Though the model does not give predictions on pt— 
dependence, some conclusions can be obtained from its 
basics (or from Gribov's Reggeon calculus in general). 
As it is already mentioned, the multi-particle production 
is described in terms of cut-pomeron diagrams. On the 
other hand, using Eq. ([T]) it is easy to see that the mean 
number of cut-pomerons increases with energy ^ s^. 
Average transverse momenta of produced particles in- 
crease with the number of chains: pt ^ \/N . So, we 
expect at LHC energies to have mean transverse momen- 
tum of produced particles higher then the one measured 
at ISR energies. Analogously, the selection of events with 
high multiplicity results in an increase of cut-pomerons 
number and to an increase oi < pt >~ ^/N . This de- 
pendence is observed experimentally by UAl and CDF 
collaborations ([15j and 16J). Only in cases of [TB] the 
values of measurements are published and in Fig. [3] we 
compare our expected dependence with these data. The 
fit result is: 



Pt = (0.62973 ±0.00084) 



(0.0717848 ± 0.00024) ^/N^. 

(5) 



In order to convert rapidity to pseudorapidity distribu- 
tions we used experimentally measured < pt > for differ- 
ent types of particles (here we assume that the sample of 
secondary charged particles consists from 7r+, and 
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Figure 1: Comparison of phenomenologically expected 
dependence of mean transfers momentum of secondary 
particles on charged particles multiplicity with experi- 
mental data from CDF 1161. 



protons, and their anti-particles). In Ref. [8] we have 
parameterized < > as a function of -y/s for charged 
pions, kaons and anti-protons, and found a successful 
fit to the data from ISR to Tevatron energies with the 
second-order polynomial function of ln(s). Fit result is 



0: 



<Pt >= 0.34 - 0.002 In s + 0.00035 In s, 

< >= 0.55 - 0.031 In s + 0.001 In^ s, 

< >= 0.65 - 0.045 In s + 0.0036 In^ s. 



(6) 



In the following for each energy and for each particle 
type we calculate pt using 

In Fig. [2] we give description of SppS and Tevatron data 
on charged particles pseudorapidity distributions in pp 
non-single diffractive (NSD) events and give a prediction 
for =14 TeV. 

Energy dependence of charged particles pseudorapid- 
ity density for SppS NSD events in the central rapid- 
ity region is compared with data in Fig. [3l In the su- 
percritical Pomeron (A > 0) theory with account of 
"non-enhanced" diagrams (without interaction between 
Pomerons), which we consider now, inclusive cross- 
sections dac/dy at very high energies and at y ~ in- 
crease with energy as (s/sq)^ (see the first article in 
Ref. [1] for more details). This means, in particular, 
that an energy dependence of inclusive spectra in the 
central rapidity region gives more reliable information 
on the value of A, than an energy dependence of crtot, 
where pomeron cuts strongly modify energy dependence 
compared to the pole diagram. An integral over ra- 
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Table 1: Predictions for LHC on mean number of charged particles in two pseudorapidity bins and in full phase- 
space . 



y/s TeV 






iV^^^'^(|r;| < 2.5) 


^7ne((-|^| ^ 2.5) 
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Figure 2: Comparison of model's prediction with data 
on charged particles pseudorapidity distribution in NSD 
events and prediction for LHC: a) description of UA5 
data [17], b) description of CDF and P238 data (circles 
and squares, respectively) [18]. 



pidity density gives charged particles mean multiplicity 
< Nch >■ In Tables [1] and [21 we list QGSM predictions on 
pseudorapidity density at mid-rapidity and mean num- 
ber of charged particles for pp interactions at various 
energies of LHC. Charged particles mean multiplicity is 
calculated for |?7| < 1 and \ri\ < 2.5 pseudorapidity bins, 
corresponding to ALICE [20] and ATLAS/CMS [21] de- 
tectors central barrel acceptances, respectively, and for 
full phase-space. 
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Figure 3: Energy dependence of charged particles pseu- 
dorapidity density at = for NSD events. Data are 
taken from [ITl [111 [19] . 



Table 2: Predictions for LHC on charged particles pseu- 
dorapidiy density at mid-rapidity. 



TeV 




(dA^^"^7d?7)„=o 


0.9 


3.5 


3. 
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4.9 


4.2 


10 


5.2 


4.5 


14 


5.5 


4.7 



The multiplicity distribution in the model is given ac- 
cording to Eq ([2]) by a sum of n cut-pomerons contribu- 
tions and it is assumed Poisson-like form for each con- 
tributor. In this scheme, the KNO scaling [22] is ap- 
proximately valid up to ISR energies {^/s < 100 GeV, 
where the Poissonian distributions for a different num- 
ber n of cut-pomerons significantly overlap) and must be 
definitely violated at higher energies [21], as confirmed 
by measurements done at SppS and Tevatron. In Fig. |4| 
we give a description of UA5 data [25 on charged par- 
ticles multiplicity distribution at -^/s = 900 GeV and 
for various pseudorapidity intervals. In Fig. [5| predic- 
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tions for ^/s = 14 TeV for the full phase-space and for 
two preudorapidity bins (corresponding to ALICE and 
ATLAS/CMS detectors central barrel acceptances) are 
given. 




Figure 4: Comparison of model predictions with UA5 
data on charged particles multiplicity distribution in 
NSD events at = 900 GeV. 
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Figure 5: Predictions for charged particles multiplicity 
distributions at LHC. 

In this formulation, the model at very high ener- 
gies (starting from ^/s ^ 10 TeV) predicts s^/ln^s 
asymptotic behavior for particles rapidity density at 
mid-rapidity and s'^/lns for mean number of secondary 
particles in the full phase-space. We stress that the 
contribution of enhanced diagrams (the interaction be- 
tween Pomerons) that we do not take into account here 
are expected to be essential at these energies. With 
the account of these diagrams {dN/dy)y=o ~ Zn^s and 
< N >~ In^s asymptotic behaviors are expected. MC 



realization of the model with account of enhanced dia- 
gram is formulated in |25j . 

In Fig. [5] we give a description of UA5 data on 
charged particles pseudorapidity distribution in single 
diffractive events at y/s = 900 GeV (dissociation of one 
of the colliding particles is considered) and give a predic- 
tion for -^8=14 TeV. For both energies the spectra are 
calculated by integrating over all masses of the diffrac- 
tive system (up to M^/s < 0.05). More comparisons 
with UA4 and UA5 data on SD events can be found in 

i- 




Figure 6: Charged particles pseudorapidity distribu- 
tions at = 900 GeV (solid line) and 14 TeV (dotted 
line) in singe diffractive dissociation events. The data 
points are from UA5. The indicated errors are statisti- 
cal and the systematical errorsareunknown. 



4 Conclusion 

In this paper we compare QGSM predictions with SppS 
and Tevatron data on charged particles pseudorapidity 
and multiplicity distributions and give predictions for 
LHC. We stress that there are no free parameters in this 
analysis and hope that this approach will give a reliable 
predictions for particle production at LHC energies. 
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